Semi-fluxons in long Josephson O-rc-junctions. 
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We investigate analytically long Josephson junctions with phase ^-discontinuity points. Such junctions are 
usually fabricated as a ramp between a superconductor like YBa2Cu307 with d-wave symmetry of the order 
parameter and an s-wave superconductor like Nb. From the top, they look like zigzags with Jt-jumps of the 
Josephson phase at the corners. These Jl-jumps, at certain conditions, lead to the formation of half-integer 
flux quanta, which we call semi-fluxons, pinned at the corners. We derive a version of sine-Gordon equation 
which describes the dynamics of the Josephson phase in such structures, and obtain an explicit formula which 
describes the shape of a semi-fluxon. Some properties of semi-fluxons are discussed. We propose a way to 
construct artificial 7t-junctions using only .s-wave superconductors. 

PACS numbers: 74.50.+r, 85.25.Cp 74.20.Rp 



I. INTRODUCTION 

Recent experiments with YBa2Cu307-Nb ramp long 
Josephson junctions (LJJ) fabricated in a zigzag geometry (if 
viewed from the top) clearly demonstrated that due to the spe- 
cific order parameter symmetry the LJJ consists of alternating 
facets of 0, 71, 0, 71... junctions!!! As a result, half-integer 
flux quanta can be spontaneously generated and trapped at 
the corners of the zigzag, because these are exactly the points 
where the order parameter of high-X c superconductor changes 
its sign due to a 90° change in the direction of the loseph- 
son contact (direction of Josephson tunneling current)!. Half- 
integer flux quanta, further called semi-fluxons (SFs), were 
also expemnentally observed in tri-crystal grain boundary 
(GB) LJJsEdQ. The presence of alternating 0- and 7C-facets 
results in a set of 7C-discontinuities of the Josephson phase at 
the corners where 0- and 7C-facets join. The possibility to fab- 
ricate such LJJs opens new perspectives for Josephson elec- 
tronics (digital circuits, fluxon devices, quantum bits, etc.p, 
as it removes certain limitations of conventional circuits, e.g., 
allows to build RSFQ-like circuits with minimum -number of 
bias resistors which means much lower dissipationEl. 

In section II we derive the version of sine-Gordon equation 
which describes the dynamics of the Josephson phase in LJJ 
with alternating and tc regions (facets). This equation de- 
scribes all possible excitations, such as fluxons, semi-fluxons, 
plasma waves, etc. In section III we obtain an explicit expres- 
sion for the semi-fluxon shape and shortly discuss the prop- 
erties of semi-fluxons. Section IV concludes this work and 
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FIG. 1 : Sketch of a fragment of a zigzag LJJ (a) and its equivalent 
circuit (b). 



only difference between a conventional LJJ model and our 
model is the. presence of phase 7C-discontinuity points along 
the junctions. 

The Kirchhof equations for currents [Fig. |](b)] and for the 
Josephson phases in an elementary loop are 
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(x+dx)-$(x) = —[® e -I L (x)L(x) 
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I L (x) +I e (x) = I L (x+dx)+I(x), 



■n(x), 



(i) 

(2) 



presents some ideas on future investigations of such junctions. 



n. DERIVATION OF THE BASIC EQUATION 

We consider a zigzag LJJ as a planar ID LJJ with a curve- 
linear coordinate x along the zigzag, as shown in Fig. |l](a). 
A magnetic field is applied perpendicular to the plane of the 
structure. Such a model can be used to describe both GB 
and ramp zigzag LJJs. For ramps, the direction of magnetic 
field can also be different, and one has to calculate the effec- 
tive field following the guidelines presented in Ref. ^| The 



where §(x) is the Josephson phase at point x of the junction, 
&e(x) is the external magnetic flux applied to the cell, L(x) is 
the inductance of the piece of the junction electrodes between 
x and x + dx, Il(x) is the current in the electrodes, i.e., through 
the inductance L(x), I e {x) is the externally applied bias cur- 
rent, and I(x,t) is the current through the Josephson junction. 
The particular expression for I(x,t) depends on the J J model 
adopted and is introduced later. 

The function TI(x) can be equal to or ±7E and shows 
whether there is a corner (discontinuity point) on the interval 
from x to x + dx. The values ±jt describe the direction of the 
jump. The function Yl(x) acts as there would be a generator 
of phase ±7C inserted in the corresponding elementary loops. 
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Imagine that we have a function Q(x) which is constant every- 
where and jumps by ±71 at each corner. This function can be 
written as 



k=\ 



(3) 



where Ok = ±1 defines the direction of the k-th jump and the 
sum is over all N c corners located at x = Xk, H(x) is the Heav- 
iside step function. 

Assuming that dx is smaller than the distance between the 
corners, the function II can be written as follows 



U(x)=Q(x + dx)-Q(x), 



(4) 



Assuming that the interval dx is infinitesimal, we can 
rewrite Eqs. ([[]) and (||) in a differential form using the fol- 
lowing expressions: 

/ = j(x)wdx, (5) 
h = j e (x)wdx; (6) 

L = ^dx; (7) 
w 

4> e = ^o(H- n)Adx = /JoH(x)Adx, (8) 

where fjQd' is the inductance of one square of the supercon- 
ducting electrodesEILaf' ss 2Xl is the effective magnetic thick- 
ness of the junctioro, n is the unit vector normal to the plane 
of the junction cell as shown in Fig. |l|(b), A w 2Xi is the effec- 
tive penetration depth of the magnetic field into the junctionlHJ, 
A-l is the London penetration depth of the superconducting 
electrode, and w is the width of the junction, e.g., for a GB 
LJJ it is equal to the film thickness. We assume that the films 
are spatially uniform so that w, d' and A are independent on x. 

Substituting Eqs. into Eqs. ([]]) and (^) we rewrite 

them in a differential form (dividing by dx — > 0): 
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HA- 



dh , . 

~T = {Je-j)W. 

dx 



(9) 
(10) 



Here and below, the subscripts t and x denote the partial 
derivatives with respect to time t and coordinate x, respec- 
tively. 

Excluding Il(x) from Eqs. ([)[) and (jici|), we get the equation 
which describes the dynamics of the Josephson phase in the 
system 

Ue -j) = ^ {"oAff* (X) - ^ [♦« - 6„(*)] } > (ID 

For the Resistively Shunted Junction model, one should sub- 
stitute 

;(*) - jcM® + p-$t +C'p-$ tt (12) 
2np 2k 



into Eq. (|TT]). Here j c , p and C' are the critical current den- 
sity, specific resistance and specific capacitance of the junc- 
tion, respectively. After this, Eq. ( |TT| ) can be rewritten in a 
form which resembles the usual sine-Gordon equatioro: 



kjtyxx - (Hp %t - sin(<|>) = co c % - y(x) 



QH x (x)+Xjd xx (x), 



(13) 



where Xj = ^/<$>q/ (2%/jq j c d') is the Josephson penetration 

depth, <Mp = \/2k j c / (OoC) is the Josephson plasma fre- 
quency, co c = 2jc j c p I <Pq is the characteristic frequency, y(x) — 
je(x)/j c is a normalized bias current density, and Q = 
2tihqAX 2 j/^q. 

For theoretical investigation of the system we introduce 
standard normalized units, i.e., we normalize the coordinate to 
the Josephson penetration depth Xj, and the time to the inverse 
plasma frequency CO^ 1 . After such simplifications, Eq. ( pj| ) 
can be rewritten as: 

<|>xt - §tt - sin((j)) = atif, - y(x) + h x (x) + d xx (x), (14) 

with the damping coefficient a = CO p /co c = 1/ a/Pc, and the 
field h normalized in the usual way as h(x) — 2H{x)/H c \, 
where H c \ = / (%^qAXj) is the first critical field (penetra- 
tion field) for a LJJ which is, in fact, equal to the field in the 
center of the fluxon. From now on all quantities are given in 
normalized units. 

In comparison with the usual perturbed sine-Gordon equa- 
tion, Eq. (JM) contains an additional term 



d xx (x) = 7C £ Gk&x (x-Xk). 

k 



(15) 



which describes the corners with ±7C phase jumps. 

To simplify the analysis it is convenient to present the phase 
(j) as a sum of two components: the magnetic one fi(x) and the 
order-parameter related one 8(x) O), i.e., 



ty(x,t)=fi(x,t) + Q(x). 



(16) 



In this case we can get rid of 8-functions in the Eq. ( 14) and 
rewrite it only for the "magnetic" component /j: 



-H„ - sin(^)cos(0) = afi t -y(x)+h x (x). 
±i 



(17) 



It is rather interesting that this is just the usual perturbed sine- 
Gordon equation, but the sign of sin(/j) changes from facet to 
facet. This means that every second facet can be considered 
as having negative critical current —1 (in normalized units) 
instead of +1. Note, that this is only valid for a current-phase 
relation with odd harmonics. In the general case odd harmon- 
ics change the sign, and even harmonics do not. This applies 
to both sine and cosine harmonics in the Fourier representa- 
tion of the current-phase relation. 

It may be easier, especially for numerical investigations, 
to use Eq. ( |l7| ) instead of Eq. (p"4|). Eq. ( ]l7| ) can be solved 
separately on each interval between corners, and all solutions 
should be joined at x — x^. On the other hand, solving Eq. ( |l7| ) 
implies dealing with 8-functions, which may be rather cum- 
bersome. 
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III. SEMI-FLUXON 

As was found experimentally!, the presence of K discon- 
tinuities of the phase, may result in the formation of semi- 
fiuxons pinned at the corners of the zigzag. Let us consider 
an infinitely long LJJ with a single corner at x = and derive 
an analytical expression which describes the shape of such a 
semi-fluxon. We start from the static version of Eq. (|l7|) with- 
out perturbation terms 



±sin/j, 



(18) 



where we have assumed that Q(x) = — nH(x), i.e., phase jumps 
from to —71 when we pass the corner at x = 0. /j~ (x) refers to 
the left half of LJJ (x < 0), while fJ + (x) refers to the right half 
(x < 0). The semi-fluxon is generated to compensate a phase 
jump at x = 0, and, far from the corner, the LJJ should not 
"know" about the jump. Therefore, we search for a solution 
of Eq. dl8[ ) which has the following boundary conditions at 
infinity 



(j)(±°°) — 0; and (j) x (±°o) = 

The same conditions for fj(x) are 

°°) = 0; = 7t; 

A»,(±°°) = 0. 



0. 



(19) 

(20) 
(21) 



We multiply both sides of Eq. ( |T8| ) by 2/jJ and rewrite it in 
the form 



= =f2(cos/j =f ) 



After integration we get 



(pj) =T2co SjU:F +C. 



(22) 



(23) 



The integration constant C can be determined from the condi- 
tions (M). Taking the limit of Eq. (^3j) at x — ► ±°° we see that 
Eq. (|23[) holds only provided C = 2. Thus 

( iU ;) 2 =2(l-cos iU -)=4sin 2 ^-; (24a) 

( ( u+) 2 =2(l+cos i u+)=4cos 2 ^. (24b) 

Let us introduce a new variable — fJ ± /2 and take the 
square root of both parts of Eq. (Q). We suppose that /jf > 
for all x, i.e., fj grows from to K always with non-negative 
derivative. This can be checked later when we will get a so- 
lution. Therefore we keep only the plus sign in front of the 
square root (the minus sign corresponds to a negative semi- 
fluxon), so we get 



\\i x = sin\|/ ; 
\|/+ = cos\|/ + . 

Integrating this equation yields 

d\\i 



sin\|/~ 

cosy 4 



= In tan — . 

2 



= ln 



1 +sin\|/ 4 

COS\|/+ 



(25a) 
(25b) 



(26a) 
(26b) 



Using the condition \|/(0) = /j(0)/2 = Jt/4, we can determine 
the value of the integration constant x* : 



x„ = In 



tan^ = In {^/l- = \ng . 



(27a) 



xi = lni±^=ln(V2+l)=lnI, (27b) 
* cos(7t/4) V / g 

where Q = tan(7t/8) = \/2- 1 w 0.404. 
Solving the Eq. ( |26| ) for \\t = fj/2 and using (^) we get 



/j (x) = 4arctan {Qe x ) . 

+ ^ a 1 - Q e ~ x 

u (x) = 4arctan 

= 7C-4arctan(^e _x ) 



(28a) 



(28b) 



The final expression for the semi-fluxon shape in terms of 
the total phase §(x) can be written in a more compact form as 



-4sign(x) arctan ( Q e 



(29) 



If we want to calculate the magnetic field we should not 
forget that (j)(x) consists of two components: (1) Q(x) phase 
jumps at the corners, and (2) the magnetic component /j(x) 
which describes the smooth variation of the phase. It is the 
derivative of fi(x) that is equal to the magnetic field at a given 
point of the junction. From Eq. (Q) the field fi x (x) is given by 
the expression 



cosh(|x| — In g ) 
The field in the center of the semi-fluxon is 
2 4 

«*(0) 



coshln^ g 



V2, 



(30) 



(31) 



and should be compared with the field in the center of a fluxon, 
which is equal to 2. 

The supercurrent density can be calculated as 

. , sinhflxl — In G ) 
sin <» = = -2sign (jt J \ »\ . (32) 
cosh (| x | — Ing) 

The functions §(x), fi(x), fi x (x) (magnetic field), and sin((j)) 
(supercurrent) are shown in Fig. ||. The very difference be- 
tween a fluxon and a semi-fluxon is that (a) the fluxon carries 
one quantum of magnetic flux, while the semi-fluxon carries 
only half of the flux quantum (therefore the name), and (b) the 
semi-fluxon has a sharp maximum which looks like a cusp. It 
would be very interesting to compare the shape of a SF ob- 
tained by scanning SQUID microscopycl with the shape given 
byEq.©. 



The Eq. ( 29 ) describes a positive semi-fluxon (PSF), i.e. the 
one containing +<t>o/2. To describe a negative semi-fluxon 
(NSF), one just have to alter the sign in front of Eq. ( p^ ) or 
to change x — > — x in Eq. (]28|). In the same time, one should 
keep in mind that the sign of 9 [of Ck in Eq. (Sb] is in no way 
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FIG. 2: Comparison of fluxon and semi-fluxon shapes, (a) The 
behavior of total phase <])(x) and of magnetic component fj(x) only, 
(b) Magnetic field profile [J x (x). (c) Supercurrent profile sin(<])) = 



related to the polarity of the SF. One can as well construct 
the PSF which sits at the point where 9 jumps up from to 
+K. In this case the total phase twist will be equal to 271, but 



physically the situation will not change. Shortly, it is only the 
sign of /v, but not [Ok in Eq. (|^)] which defines the polarity 
oftheSF. 

SFs are very similar to fluxons when they interact with each 
other: semi-fluxons and (semi-)fluxons of the same polarity 
repel each other, while the ones of opposite polarity attract 
themselves. This can be easily shown writing the potential 
energy as a function of the distance between them. 



IV. CONCLUSIONS 

We have derived the perturbed sine-Gordon equation which 
describes the dynamics of the Josephson phase in a LJJ con- 
taining phase ^-discontinuities, which correspond to the cor- 
ners of the Nb-YBCO zigzag LJJ. Using the derived basic 
Eq. ( p"4] ) we have obtained the shape of a semi-fluxon — the 
new type of object which appears due to the phase jumps. 
Our results allow to investigate the interaction between semi- 
fluxons as well as between fluxons and semi-fluxons. This is 
also a starting point for the numerical simulation of various 
aspects of fluxon and semi-fluxons dynamics. 

An interesting consequence of Eq. (jl^) is that the terms j(x) 
and Q xx (x) play a similar role in the equation. This means that 
one can substitute one by the other. For example, one may 
wish to create an artificial ^-discontinuity point using only 
superconductors with s-wave order parameter symmetry with 
an injector and receptor of current of the size Ax separated 
by a minimum distance. Passing the current equal to 4-n/Ax 2 
from the injector to the receptor, one emulates the effect of 
Q xx . Of course, Ax must be much smaller than any character- 
istic length, e.g., Ax <g; Xj. Although this emulation may look 
not ideal, one should keep in mind that in a real zigzag junc- 
tion, the size of the corner is also finite and is defined by the 
lithographic accuracy ~ 1/jm. 
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